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Abstract 

In this paper we deal with the long time existence for the Cauchy problem associated to some 
asymptotic models for long wave, small amplitude gravity surface water waves. We generalize 
some of the results that can be found in the literature devoted to the study of Boussinesq systems 
by implementing an energy method on spectrally localized equations. In particular, we obtain 
better results in terms of the regularity level required to solve the initial value problem on large 
time scales and we do not make use of the positive depth assumption. 

Keywords Boussinesq systems; Long time existence; 


1 Introduction 

1.1 The abed systems 

The following abed Boussinesq systems were introduced in [5] as asymptotic models for studying long 
wave, small amplitude gravity surface water waves: 

J {I — ebA) dtT] + divV + aediv AV+ £ div (rjV) = 0, 

\ (/— edA) + V 77 + ceVA ?7 + eiy |P|^ = 0. 

In system (HH) e is a small parameter while 

/ r] = r]{t,x) gR, 

\ V = V{t,x) G K", 

with {t, x) G [0, 00 ) X R" are approximations of the free surface of the water and of the fluid’s velocity 
respectively. As it will soon be clearer, we mention that the only values of n for which (II.1|) is 
physically relevant are n = 1,2. The above family of systems is derived from the classical mathematical 
formulation of the water waves problem: considering a layer of incompressible, irrotational, perfect 
fluid flowing through a canal with flat bottom represented by the plane: 

{ix,y,z) : z = -h}, 
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where h > 0 and assuming that the free surface resulting from an initial perturbation of the steady 
state can be described as being the graph of a function 77 over the flat bottom, the water waves problem 
is governed by the following system of equations: 


A(j) + d^fj) = 0 in 
= 0 at 

diTj + V(j)Vri — dz4> = 0 at 

+ 5 Z = 0 at 


-h< z< r]{x,y,t) , 
z = —h, 
z = v{x,y,t) , 
z = f]ix,y,t) 


(WW) 


where (j) stands for the fluid’s velocity potential and g is the acceleration of gravity. The operators A 
and V are taken with respect to {x,y). Of course, in many applications the above system of equations 
raises a significant number of problems both theoretically and numerically. This is the reason why 
an important number of approximate models have been established, each of them dealing with some 
particular physical regimes. The abed systems of equations deals with the so called Boussinesq regime 
which we will explain now. We consider the following quantities: A = ma-Xx^y^t \'n\ the maximum 
amplitude encountered in the wave motion, I the smallest wavelength for which the flow has significant 
energy and cq = the kinematic wave velocity. The Boussinesq regime is characterized by the 
parameters: 



( 1 . 2 ) 


which are supposed to obey the following relations: 


a 1, /3 1 and S' ~ 1 . 

Supposing actually that S = 1 and choosing e = a = P, the systems dm are derived back in [5] . The 
unknown functions ( 77 , F) in ( 11 . 11 ) represent the deviation of the free surface from the rest state while 
y is an O (e^) approximation of the velocity taken at a certain depth. Actually, the zeros on the 
right hand side of dm represent the O (e^) terms neglected in establishing (HD- The parameters 
a, &, c, d are also restricted by: 

a-\-h-\-c-\-d = —. (1.3) 

Asymptotic models taking into account different topographies of the bottom were also derived, see for 
instance [3], Section 2, for bottoms given by the surface: 

{(x, y,z) : z = -h + eS{x, y)}, 

where S is smooth enough or [5] for slowly variating bottoms i.e. the function S = S {t, x, y) depends 
also on time. A systematic study of asymptotic models for the water waves problem along with their 
rigorous justiheation can be found for instance in nj. 

The study of the local well-posedness of the abed systems is the subject under investigation in 
several papers see for instance 0 . for space dimension n = 1 or |^, [^ (the BBM-BBM case 
6 = d > 0, a = c = 0), [10] (the KdV-KdV case 6 = d = 0, a = c > 0) for dimension n = 2. In [3] it 
is shown that the linearized equation near the null solution of dm is well posed in two generic cases, 
namely: 


a<0, c<0, 6>0, d>0 (1.4) 

or a = c > 0 and 6 > 0 , d > 0 . (1-5) 

It was proved in |7] (see also m) that the error estimate between the solution of dm and the water 
wave system is cumulating in time like O (s^t). Thus, solutions of (11.11) that exist on time intervals of 
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order O (i) are good approximation for (IWWI) as the error remains of small order i.e. O (e). Actually, 
due to the previous mentioned error estimate, on time scales larger then O (^) the solution {r],V) 
stops being relevant as an approximation for the original system. The question of long O (i)-time 
existence of solutions of dm did not receive a satisfactory answer until in |12) where, the case (d 
was treated and long time existence for the Cauchy problem was systematically proved, provided that 
the initial data lies in some Sobolev spaces. The difficulty comes from the lack of symmetry (the 
erj div V term from the first equation of the abed-system) of (11.111 . Because of the dispersive operators 
—ebAdt + odiv A, —edAdt + eVA, classical symmetrizing techniques for hyperbolic systems of PDE’s 
fail to be successful. 

Global existence is known to hold true for (ED in dimension n = 1 for some particular cases: the 
case 

a = b = c = 0, d> 0, 

studied by Amick [T] and Schonbek m and in the case 


b = d> 0, a < 0, c < 0, 

assuming some smallness condition on the initial data, see [5]. In both cases, it is assumed that: 

inf {l + e? 7 o(a;)} > 0 

a condition that makes perfect sense from a physical point of view as 1 + £77 (t, x) represents the total 
height of the water above the bottom in x at time t . The latter case uses the Hamiltonian structure 
of the system, namely, when b = d we have: 

J + (1 + erj) V'^ — ec [dxT])^ — ea dx = 0. (1.6) 



1.2 The main result 

The aim of this paper is to generalize most of the results presented in [12], namely, we address the 
long time existence issue for the general abed systems. More precisely we wish to construct solutions 
for (II. ip for which the time of existence is bounded from below by a O (i)-order quantity. The key 
ingredient is that we establish energy-type estimations for spectrally localized equations and by doing 
so we require a lower regularity index than the one used in HH. Also, we avoid using the non-cavitation 
condition 

1 + svo (x) > a > 0, (1.7) 

imposed on the initial datum rjQ or the curl free condition on the initial data Vq used in [Sj. As opposed 
to the method used in |12| , ours permits us to treat in a unified manner most of the cases corresponding 
to the values of the o, b, c, d parameters. In order to carry on our approach we need some restrictions 
on the a, b, c, d parameters. More precisely, we will consider the case (ED such that b-\-d > 0 excluding 
the following two cases: 

f a = d = 0, c < 0, & > 0. . . 

|^a = & = 0,c<0,(i>0. 

In Section [5. II we put forward the basic ingredients in order to obtain long time existence for the former 
two cases. In view of (|1.3L & -|- d > 0 is not restrictive as far as a, c are less or equal to 0. 

First, we will slightly change the form of ED noticing that the divergence free part of V remains 
constant during time. Indeed, formally, if (77, H) is a solution of E with initial data 

77|t=o = Vo I %=o = ^^ 0 , 
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then 


dtV =—V {I — edA) ^ [ r] + csAr] + e-\V\ 


and consequently we have that: 
where 


dfFV = 0 , 


is the Leray projector over divergence free vector fields. Thus, we get that: 


TV = Pt4, =‘ IV. 


(1.9) 


Of course, we have that 

div W = 0. 

By setting 

r Fo = Vj, + IT , 

I v=v+w 

we infer that the system verified by the couple (77, T) is the following: 

{ (/ — sbA) dtT] + div V + ae div AT + elTV?7 + e div (77T) = 0, 

(/ - sdA) dtV + V?7 + ceVAri + e^V |1T|^ + eVTTT + eVTlT + eiV |T|^ = 0, (1.10) 

77|t=o = Vo , y\t=o = Vo- 

Also because of m we get that curlT = 0 at any time meaning that for all I, k G 1, n we have: 

diV>^ = dkVK ( 1 . 11 ) 


The advantage of working with system (jl.l0|l instead of (II.Ill is two-folded. On one side certain 
commutators involving V behave better if its curl is 0 and on the other side, for some values of the 
a, b, c, d parameters e.g. a = c = d = 0, 6>0we need less regularity on the initial data Vq. 

Let us establish some notations. The spaces (R") with p G [l,oo] will denote the classical 
Lebesgue spaces. Given s G R we will consider the following set of indices: 

{ Si = s -(- sgn (b) — sgn (c), 

S2 = s -(- sgn (d) — sgn (a), (1-12) 

S3 = s -I- 1 - sgn (a ), 

where the sign function sgn is given by: 

( 1 if a; > 0, 

sgn (a;) = < 0 if a; = 0, 

[ —1 if a; < 0, 

and a, b, c, d are chosen as in dn. We will denote by iJ® (R”) the classical Sobolev space of regularity 
index s with the norm 

\\v\\ls= J (l + |?|^)'|^(Ol^de (1.13) 


4 





New long time existence results for a class of Boussinesq-type systems 


For any vector, matrix or 3-tensor field with components in (R") we denote its Sobolev norm by 
the square root of the sum of the squares of the Sobolev norms of its components. For any pair 
{rj, V) G (R") X (R”))" we will use the notation 


\\{V,V)\\1 = WvWh. +e{b-c) ||V77||^. + e^{-c)b\\V^r]\\]^, + 

+ \\V\\l. + e{d - a) \\yV\\l. + e^{-a)d\\V^V\\l^ , 

where = (af,ry). . and = (dlV'^). . Clearly, (R’^) x (R”))” endowed with IK-, ■)L 

is a Banach space which is continuously imbedded in (R") x (R"))". 

Our approach is based on an energy method applied to spectrally localized equations. We first 
derive a priori estimates and we establish local existence and uniqueness of solutions for the general 
abed system. Before we state the main result let us formalize the notion of long time existence of 
solutions for (11.101) in the next definitions. 

Definition 1.1. Let T > 0 a positive real number. We will say that T is bounded from below by a 
O (i)-order quantity if there exists another positive real number C, independent of e such that: 


T > 


C 

s 


Let us consider a Banach space {X, IHIx) which is continuously imbedded in (R") x (R"))". 

Definition 1.2. Let us consider W G (R™))". We say that we can establish long time existence 
and uniqueness of solutions for the equation (I1.10|) in X if for any (ryo, Vq) G X there exists a positive 
time T > 0, an unique solutio'^ (iq, V) € C ([0, T ], X) of (11.101) and a function F : (0, -|-cx)) —> (0, -|-c») 
independent of e such that: 

F{\\{iqo,Vo)\\^) 

~ e 

Remark 1.1. Of course, the function F appearing in the preceding definition is allowed to depend on 
a,b,c,d,W and on the dimension n. 


We are now in the position of stating our long time existence result: 

Theorem 1. Let a, b, c, d as in m excluding the two cases dUl), & -|- d > 0. Let us consider s such 
that s > ^ -|- 1 with n > 1. Let us also consider Si, $2 and S 3 defined by (11.121) and W G 
Then, we can establish long time existence and uniqueness of solutions for the equation (11.101) in 
F[‘^^x[Fl^^y^. Moreover, if we denote by T {rjo, Vq) , the maximal time of existence then there exists some 
T G [0, T (ryo, Vb)) which is bounded from below by an O (J-)-order quantity and a function G : R —>■ R 
such that for all t G [0, T] we have: 


||(ry,F)||, <G(||(7yo,Co)|lJ, 
where G may depend on a, b, c, d, s, n but not on e. 

Theorem [T] is the consequence of a more general result that we obtain later in this paper. In fact, 
our method enables us (without any extra effort) to establish long time existence and uniqueness of 
solutions in the more general Besov spaces, thus achieving the critical regularity s = f-1-1, see Theorem 

The rest of the paper is organized as follows. In Section[5]we establish all the basic energy estimates 
that we will need in order to prove Theorem [1] In Section |3] we prove that (ll.IOl) admits an unique 

^The solution should be understood in the tempered distribution sense. See Definition l3.ll 
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solution and we establish an explosion criterion. The method used to construct the solution assures 
an existence time that is bounded below by a quantity of order O • Finally in Section[4]we prove 
Theorem [T] namely we show that the solution of (II.IOII constructed in Section |3] persists on a time of 
order O (i). The proof will be a by-product of some rehned energy estimates that we prove in Section 
[5] and the explosion criteria established in Section |31 In Section 15.11 we discuss about the cases (HH). 
We end the paper with Section fS. 21 where we discuss the possibility of applying our method to the abed 
systems in the case of a general bottom topography derived in [9], Section 2. 

1.3 Notations 

Because our proof makes use of elementary tensor calculus let us introduce some basic notations. For 
any vector held U : K" —>■ M" we denote by VC/ : K" —>■ At„(K.) and by V*C/ : K" —>■ At„(K.) the n x n 
matrices dehned by: 


(VC/)y = 5.C/^ 

{V^U)^^=d,U\ 

In the same manner we dehne V^C/ : R” —>■ R" x K" x K” as: 

We will suppose that all vectors appearing are column vectors and thus the (classical) product between 
a matrix held A and a vector held U will be the vectoJl: 

{Auy = A,jW. 

We will often write the contraction operation between V^C/ and a vector held V by 

(V^C/ : V)^. =dlU'^V^ 

If C/, V : R” —>■ R" are two vector helds and A, B ■. R" —7\d„{R) two matrix helds we denote: 

UV = WV\ A: B = AijBij, 

{U,V)^. = I WV\ {A,B)^, = I A,B,, 

\\U\\l, = {U,U)^,, \\A\\l, = {A,A)^, 

\\V^U\\l^ = JvU :VU = j 

Also, the tensorial product of two vector helds U, V is dehned as the matrix held U ®V : R" —>■ Ad„(R) 
given by: 

(C/0 = C/WT 

We have the following derivation rule: if u is a scalar held, U, V are vector helds and A : R" —AIn(®.) 
then: 


V div (uU) = V^uU + Vm div U + VUVu A uV div U 
V{UV) =VUV + VVU 

^From now on we will use the Einstein summation convention over repeted indices. 
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ViAV) = : F +VyA* 

If we suppose that curl 1^ = 0 then the following integration by parts identity holds true: 

{VV : U,VV)^2 = j (VV ■.U):VV = j dfjV^U'^d.V^ = j 

Let C be the annulus G R" : 3/4 < |^| < 8/3}. Let us choose two radial functions y G 
P(i3(0,4/3)) and ip G ViC) valued in the interval [0,1] and such that: 

V^GR", y(C) + ^:^(2-^'0 = l. 

j>0 


Let us denote hy h = J- and h = T ^y. For all u G S', the nonhomogeneous dyadic blocks are 
defined as follows: 


AjU = 0 if j < —2, 

A-iu = xiD) u = h* u, (1-15) 

AjU = ip (2~^D) u = 2^"^ j h (2‘^y) u{x — y) dy if j > 0. 

The high frequency cut-off operator Sj is defined by 

SjU — AkU. 

k<j-l 


Let us define now the nonhomogeneous Besov spaces. 


Definition 1.3. Let s G R, r G [l,oo]. 
u G S' such that: 


ll^ll 


The Besov space i3| ^ 





is the set of tempered distributions 
< oo. 


Let us mention that LT^ (R") = 2 (M") and that we have the following continuous embedding 


Bh Bl^ ^ H\ 

for all s' < s. Some basic properties about Besov spaces can be found in the Appendix. For more 
details and full proofs we refer to [H]. 


Let us consider e < 1 and s > 0, r G [1, 00 ]. For all ( 77 , V) G (1^") x (R"))” we introduce 
the following quantities: 

U] = U] ( 7 ?, y) = y 7?| + e (& - c) |Vr;,f + e\-c)h {V\ : V\) (1.16) 

+ Jvf + eid-a) {VV, : VV,) + e\-a)d {V^V, : V^V,) , 


and 
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U^ = U!irj,V) = M% +eib-c)\\Vr,\\l. + s^-c)b\\V^r,\ 


(1.17) 


+ e{d-a)\\VV\\ls +e'^i-a)d\\V^V\ 


where {r]j, Vj) := (A^ry, AjV) are the frequency-localized dyadic blocks defined by relation (11.151) . It is 
easy to check that Ug {rj, V) is a norm on the space (K”) x (K")) . Also, it transpires that: 




Ug. 


Some times we will also use the notation 

mV)\\^ = Ug{rj,V), 

for all {r,,V) e B 2 .V (R”) x (^ 2 (r(®"))" ■ We observe that (B 2 .V (»") x (^ 2 ),. (M”))" , IK-, OIL) is a 
Banach space. In order to ease the notations we will rather write i?| ^ instead of i?| ^ (R"). Another 
quantity that will play an important role in the following analysis is: 

N] = N] {r],V) = J il+e\\r]\\^^)Tj] +e{b-c) (1 + e llTylL-) 

+ j £^(.-c)bil+e\\'n\\^^ ){V^r]j -.V^rjj) 

+ J {l+sr] + e\\r]\\^^)Vj^+e{d-a + der] + d£\\r]\\^^){VVj :VVj) 


which satisfies: 
Denoting by 

we obviously have that 


2 Energy-type identities 

We begin by localizing equation (11.101) in Fourier space thus obtaining that: 

(I — ebA) dtr]j + div Vj + ae div AVj + eWVijj -b eWrjj + ery div Vj = eRij 
(I - sdA) dtVj + Vrjj + ceVAry^ -b eVy IF + eVVjV = £R 2 j 
Vj\t=Q = Ajijo , yj\t=o = AjVb 

where the remainder terms are given bill: 

r Ri, = [W A,] V?y + [F, A,] Vry + [ry, A,] div F 
j i?2, = [IF, AO VF + [F, AO VF - iVA, (|1FL) - A, (VIFF). 

^From now on, we agree that if A, B are two operator then the operator [A, B] is given by: 

[A,B] = AB-BA. 


0 : WO , 

(1.18) 

,=o)^C/,(ry,F) 

(1.19) 


(1.20) 

^)^Ug{v,V). 

(1.21) 


( 2 . 1 ) 


( 2 . 2 ) 
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Let us establish our first useful identity. We multiply the hrst equation in (j2.ip by rjj and the second 
one with Vj and by adding them up and integrating in space we get thal0: 

{\\Vj\\l 2 + ebWVrjjWl^ + \\Vj\\l 2 + ed || VV,- J Vj div AV^ + ce j VjV Arjj (2.3) 

+£ J ririj div = I y div V (rj^ + Vf) + e j RijTJj +e J R 2 jVj. 

Let us denote by Ti the right hand side of the above identity: 

Ti = I y div 1/ (r?f + y/) + e y R,,r,, + e J R^.V,. (2.4) 

Next, we wish to derive similar identities involving the quantities Vr]j,V^rjj and VV,, V^V^. In order 
to do so, let us observe that applying V to the hrst equation in (11.101) gives us: 

(/ — ehA) dtVr] + V div V + aeV div AV + eVWVr] + eV^-qW + sV div (viV) = 0 

and that by applying Aj we end up with: 

(/ — ebA) dtVqj + V div Vj + aeV div AVj + eV^qjW + eV^qjV + eqV div Vj (2.5) 

+edivVjVq + eWjVq = —eAj {VWVq) + eR^j, 

where 


i? 3 j- = [W Aj\ + [1/, Aj] V^q + [q, A^-] V div ^ + [V 77 , A^-] div V + [Vq, Aj] VV 
We multiply (12.5p with —ceVqj and by integration we gef@: 


—ce 


dt (^J (\Vqjf +eby^qj : ~ J divVjVqj 


— ace^ / V div AVjVTjj 


( 2 . 6 ) 


—ce^ / qV div VjVqj = T 2 


with T 2 given by 

T2 = ce^ y (V^qjW)yqj + C£^ y (y^qjV)yqj 

/div V, V,V,, + /(Vld V,)Vd, + ce‘I A, iyWWr,) V,, - 

We proceed similarly with the second equation in (11.101) and we obtain: 


(/ - edA) dtVV + V^q + ecV^Aq + \Wf + £ (V^W : V + VV"V*W) 

^Observe that here we use the fact that curl V = 0. Indeed, under we have 

/ ^VjVVj = \jy^\V\^ = -lj divV |y,f. 

^Here, we use the fact that divVF = 0. We get that: 



(2.7) 
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+e (VV : W + VWV*V) + e (VV : V + VVV*V) = 0. 

We localize the last equation and we get that: 

(/ - edA) dtVVj + + ecV^Ar^j + |W|^ + eA^ {V^W : V) (2.8) 

+eAj (VyV*W) + eVVj : W + eA^ (VWVV) + : V + eVVjV^V = sRaj, 


where: 

We contract (12. 8 |) with —aeWj and by integration we get that: 


Raj = [W, Aj] VV + [y, Aj] VV + [VV, Aj]VV. 


—ae 


dt (^J VVj : VVj + edV^Vj : - ae j \7^T]j : VV, - ace^ J V^Arj^ : \/Vj = T 3 (2.9) 


with T 3 given by: 

T 3 = ae^ J (VVj : W) : VVj + ae^ J {V% : V) : VVj + ae^ J VFj VV : 

+ae^ / + Aj : V) + A^ (VW^W)^ : WVj 

Let us add up identities (12.61) and (12.91) to get thal[^: 

—ae J y'^rjj : Wj — ce J V div VjVrjj — ce^ J 77 V div VjVrjj = T 2 + T 3 . 

Finally add up (12.31) to the last identity in order to obtair0: 

\dt (yj Vj +e{b-c) \yr]jf + e^{-c)b {V^rjj : V^rjj)^ 

+ ^^‘(/ \VA^+^(d-a) (VF, : VV,)+e^(-a)d(vV, : VV,)) 

-\-e J rjrjj div Vj — ce^ J rjV div VjVrij = Ti + T 2 + T^. 


VF - ae^ / i? 4 ,- : VF 


( 2 . 10 ) 


( 2 . 11 ) 


2.1 Refined energy-type identities 

As we have already seen in (II. 6 |) , the abed system possesses a formally conserved energy. The types of 
estimates that we establish in this section, resemble very much to this conserved energy and will be 
the equivalent of the ones obtained in [HI pages 617 and 625, with Ajij and AjV instead of 77 and V. 
Of course this is the key ingredient for obtaining long time existence results that allow initial data to 
lie in larger spaces. 

® Observe that by integration by parts we get 

—ac£^ J V div hVj : Vr}j — ace^ J V^A7]j : Wj = 0. 

^Observe that the terms ae f rjj div AV^' — ae f : VVj and ce f VjVArjj — ce f V div Vj : Vrjj are both 0. 
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Having established identity (j2.11|) we observe that the terms e f rjrjj div Vj and — ce^ / 77 V div VjVrjj 
prevent us from directly applying a Gronwall type argument and establishing long time existence. In 
order to repair this inconvenience let us multiply the second equation of (EH) with erjVj. We thus get: 


(e {I - dA) dtVj, rjVj)+ e j rjVrjjVj + ce^ J t]VA r]jVj 
= J (yv,W) I {VV,V) (vV,) + J 

Observing that the first term writes: 

{e{I - dA)dtVj,r]Vj )^2 ^ (^J sri\Vj\^ + e^driWj ■-dtv\Vj\^ 

J dtvVVj ■.VVj+e^d{dtVVj,yTj(^Vj)^2 


we write that 


iSt (^J er] + e^drjVVj : VV^ + ej r]\/r,jVj + ce^ J riVArj^Vj = T4 

with T 4 given by: 

T4 = -e^ I (VH, W) vV, -e^ J (VV^V) tjVj + j vR2jVj + 

+ 1 J 9tV J dtV^V, : VH, - e^d{dtVV„Vr] ® ■ 

We add up (12.111) and (12.121) in order to obtain: 

^dt (^Jri^ + s{b- c)\Vrj/ + e\-c)b 
+ ^dt (^J (1 + er]) IHjf +e{d-a + derj) (VV, : VV,) + s^{-a)d (V^V, : 


= To + Ti + Tz + Ta + T 4 , 


( 2 . 12 ) 


(2.13) 


(2.14) 


with 


To — e 


VrjrjjVj — ce^ 


VrjArjjVj + ce^ 


’S/vjS/rij div Vj. 


2.2 Estimates for the Tj’s 

Having established the energy identity (j2.14|l . we proceed by conveniently bounding the RHS term. 
We want to obtain a bound of the form: 

To + Ti + T2 + T3 + T4 < £P(Uj, Us) 

with P (x, y) some polynomial function with coefficients not depending on e. We suppose that 

n n 

s> — + lors= — + 1 and r = 1. 

Moreover, C > 0 will denote a generic positive constant depending only on the dimension n and on s. 
All the estimates established here are valid for the a, b, c, d parameters satisfying (|1.4I) with 5 + d > 0. 
The case: 

a = d = 0, b>0, c<0 (2.15) 
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needs more attention and we will investigate it in Section 12.41 thus for the moment we do all the 
computation supposing that case (12.151) does nor occur. 

We claim that we can bound Tq + Ti + T 2 + T 3 + T 4 in such a way that the following estimate holds 
true: 

\dt {^j Vj +e{b-c)\Vri/ + e^{-c)b 

+ ^dt (^J (1 + er]) y/ + e (rf _ a + derj) {VVj : VVj) + e^{-a)d {V^Vj : 

< eCUj {Uj {Us + HUs + U^) + + Us) {H^ + HUs + U^)} , (2.16) 


where C > 0 depends only on a,b,c,d,n and s but not on e, {cj {t))j is a sequence with {‘2^'^Cj (t))j G 
having norm 1 and finally H is defined by: 


H = llVbllsi,, + II - sgn (a) || 


Let us detail this below. Regarding Tq, we begin by observing that: 

e I VmjV, < ellVryll^^ ||ry,||^, ||R,||^, < eUps- 

Next, if a = d = 0 then 6 > 0 and we get thalU: 


— J ArjVrjjVj + ce^ j \7g\7rjj div Vj 

= —ce^ J Ar]\7r]jVj — ce^ J V^gVgjVj — cs^ J V^rjjVriVj 

< -ce^ (||Ar?||^,. ||Vr?,||^. ||R,||^. + ||VM|^„ ||V77,|1^. ||y,|l^. 

+ \\^v\\L^PW\m\p 


< eC max 




If at least one of a, d is not 0 then we write: 


We choose 


— ce^ J VrjArjjVj + ce^ J VrjVrjj div Vj 
= ce^ J iy’^riVj) Vr]j + ce^ J (VV, V 77 ) \7r]j + +ce^ J Vr]Vr]j div Vj 


< -ce^ {\Pvh^ ||V7?,||^. + llVTyll^,. ||Vy,||^. ||V7?, 


'JWL^ 


W^vWl- l|V?7,IU.||divl/,|l^,) 


< sC max 


—c —c 




pi _ 
^abcd 


ifd-a = 0 , 

maxfjE^,^) ifd-a> 0 . 


(2.17) 


®From now on, we adopt the convention ^ = 0. 
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Consequently, we have: 


To < 


(2.18) 

Next, let us analyze Ti. According to Proposition 16.71 and Proposition 16.21 we get that: 


< Cc] it) (llVTPII^.-i + llVPlI^.-i ||P|U. J , 

\R2Al- < Cc^j it) ((iIvm^IIb^- + IIvp||b.-i) ||p||b|_^ + {\\w\\bi^ + II^^IIb.v) > 


with (2l®c* (t)) € i = 1,2, with norm 1. In order to avoid mentioning each time, from now on, 

for all natural number, i G N, (c* (t)) . will be a sequence such that (2l®c* (t)) . G £''(Z) with norm 


1. Recall that H stands for the following quantity: 


H=\\W\\s, +IIVIPII 5 , -sgnia)V^\\y^W\\ 


(2.19) 


and rewrite the previous inequalities (of course the constant C = C (n, s) changes whenever necessary) 
as: 


||Ri,|li 2 <C'c](t)t/,(i7 + I/,), 

11^2,11^2 < Ccf (t) {H^ + Us [Us + H)) . 

We get that 

Ti = I / div P ( 77 I + P/) +e J Rijijj +e j R 2 jVj 

< e ||divP||^„„ (||77,||2, + ||P^.||2^) + eCcf (t) U, {H^ + Us {Us + H)) 

< Ce {Ups + cf {t) Uj {H^ + HUs + P")). 

We turn our attention towards T 2 '. 

T 2 = ce^ J {V^T]jW)VT]j + ce^ J (V^ryjP)Vryj + ce^ J div VjVrjVrjj + ce^ J VV) VryVTy^ 

According to Proposition 16.21 and Proposition 16 . 71 we get that 


( 2 . 20 ) 

( 2 . 21 ) 


( 2 . 22 ) 


||A,(VIPVr;)||^, + ||i?3,IL2 

< Ccf (t) (IIVIPII^.^^ llVrylls.v + llVW^IIsrP 

+ llVPIls-^ l|Vr?||B._^ + \\Vrj\\^.-i IIVPII^.^^ + ||V\|| 

< Cc| (f) (||VIP||5,_^ llVrylls.v + H^Ubi,, H^^Hb.v 

+ IMb|,J|VP|Ib,v + l|V77||s,^ J|P|Ib|J • 


2v) 


(2.23) 


We observe that it is here that we need the restriction (12.151) on the parameters. Indeed if a = d = 0 
and c < 0, then we only have P G ^ hence Us cannot control || VP||^s . Let us observe that: 


I{V\W)yr,, = -iy |V7?,y divIP = 0 
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and 


JiV\V)Vii, = - iy divl/. 

We infer that: 

T, = / div F, /(vr, V,) V, +ce^j A, (Vld'V,) V,, - / i!,, V.,, 

< :^s^[\\diYV\\^^ l|V7?,||i. +2||V77||i„ ||V7?,||^. ||VF,||^. 

+ Cel (t) ||V77,||^. (llVH^IIsi,, I|V^IIb,v + I|V^I1b,v 

+ hllBj Jl w|Ib,v + liv^lls,v II^IIbi J} 

{ups + c| (<) . 


< eC max 


—c —c 
b — c d — a 


We let 


Cabed = max 


—c —c 


(2.24) 


(2.25) 


b — c d — a 

Let us estimate T 3 . As above, owing to Proposition 16.71 we may bound i? 4 j in the following manner: 

-ae^ \\R 4 j\\t 2 < -ae^Cc® (t) 


B. 


:-I|VP||b, +||VP||5.-d||VP||j 


vvw 


J 


< -ae^CcI it) (||fP|ls,v 

< Cei it) UsiH + Us). 

\/d — a 


(2.26) 


Also, we can write due to Proposition 16.21 

Q ||a,v 2 + II A, (VWVV)II^, + II A, (V^w : P)||^, + ||A, (VPV‘W)||^,) 

<-aeW^{t)(\pw\\^^ ||W||b| + ||VW|||. + ||VW||b. HVPII^. 

\ii "^ 2 ,r 2,r 2,r 2,r 2,r 


—as 




< C £2 max ( —a, 




—a 
•yj d — I 


-||VfP|lB|J|VP||B, 
c| (t) (7J2 + . 


Then, using the integration by parts identity (jl.1411 we get that: 

T 3 = ae^ J {Vp : W) : VVj + ae^ J {Vp : V) : VVj 

+ae^ J VVp*V : VVj - ae^ J p : VVj 

+ae^ J (^^Aj |W|^ + A^ (VWVV) + A^ {V^W : V) + A^ {VVV*W)^ : VVj 


< -as^ ( - IldivPIl^^ + II W||^^ live, 11^. ) + 


+Ce max 


/ —a —a 


VVd — a ’d — a 


(c| (t) e-iups [H + Us) + c| (t) £-5Cj (iJ^ + 
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< eCmax ( , iUps + c] (t) Uj (H^ + HUs + U^)) (2.27) 

Vvrf — a a — aJ 

and as before, let us denote by 

Finally let us turn our attention towards r 4 . Let us write: 

T 4 = (j '4 + 7 ^ 4 , 

with 

= dtV IVjl"" + - e^d (9tVy„ V 77 ® Vj)^, 

and observe that by integration by parts and ( 12 . 211 ) 

G 4 = -e^ J (VVjW) rjVj - J (VVjV) rjVj + j r]R2jVj 

< y/ ((div(r;lF) + div(r?y))|y,f) +G£%2(t)||y,||^.||ry||^^ ^h^ + hUs + U^) 

< Ce {Uf {HUs + U^) + c® (t) UjUs {H^ + HUs + U^)) ■ (2.29) 

Let us now analyze the term B 4 . We begin with 

I J dtv\Vjf <e\\dtr]\\^^ \\Vj\\l2 < ^Uf\\dtr]\\^^ 

We have that: 



(/ — he A) ^ [(/ + aeA.) div V + e div {r](W + F))] 

^00 

< 

(/ — be A) ^ (/ + aeA) div V 

+ ^ 

L°° 

(/ - beA)-^ 

div ( 77 ( 14 ^ + y)) 

< 

(/ — &eA)~^ (/ + aeA) div V 

- +e\\Wr]\\Bs +e\\r]V\\ss 

-°2,1 


If 6 > 0 or a = 6 = 0, then because the operator (/ — be A) ^ (/ + aeA) maps to with norm 
independent of e, we get that: 

\\dtv\\L°° ^ Gmax ^1, {Us + HUs + Uj) 


o' 

II 

)— I 

a < 0 then d > 0 and we see 

that 

we have 





— ae divAy d < 

^2,1 

—ae 

||W|| <- 

■ae 

’"II 05 ,. s ]fi"- 

(2.31) 

We set 



(1,^) if5>0 





4 f 

max 

or a = 

CD 

II 

-0 

(2.32) 



max 

(i-V?) if*' 

= 0 and 

a < 0. 
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Thus, we get that 

I j dtv ([/. + HUs + U^) . 

In a similar fashion we can treat the second term of B 4 thus obtaining 


dtvVV, : VV, < (C7« + HU, + U^) , 


and we set 

(^5 _ ^ f^4 

^abcd ^ abed' 

Finally, the last term of is estimated as follows: 


-e^d{dtVV,,Vr^®V,)^. < e^d\\Vr^\\^^ \\dtVV,\\^. \\V,\\^. , 


and using 

dtVVj = - (/ - edA)”^ V (Vr^j + ceVATy^ + eVy^IF + eVVjV - ei?2j) 
we observe that 


(2.33) 


(2.34) 

(2.35) 


ed\\dtVV,\\^. <yrd{U\\^.-ce\\Vr^j\\^.+e\\W\\j^^ \\VV,\\^. 

+eniL=o ||VV,|l^.+£p2,||^.) 

< {Uj {l + H + U,) + c] {t) {H^ + HUs + U^)) , 


where 

C.V, = m» (^v^, \/^) , 

thus, we conclude that 

-e^d (9tVy„ Vry ® < eCCl,,,U, {U, [U, + HU, + U^) + 

+c) {t) Us {H^ + HUs + U^s)) ■ 

Combining estimates (|2.33l) . (I2.34|) . (12.371) we obtain: 


(2.36) 


(2.37) 


Bi < eCCabedUj {Uj(Us + HUs + U^) + (t) U, {H^ + HU, + U^)}, (2.38) 

where 

Cabed = (2.39) 

i=l,6 

is the maximum of all constants depending on a,b,c,d that appear in relations (12.171) . (I2.25F (12.281) . 

(lO^ . (1^ . (1^ . 

Thus, supposing that we are not in the case: 


a = d = 0 and c < 0, & > 0, 
we are able to successfully bound the Ti’s. 

Finally, after adding up the estimations (|2.18|) . (I2.22F (12.241) . (I2.27F (|2.29|) . (I2.38P we obtain that 
there exists a positive constant C > 0 depending only on n and s but not on e such that 

\dt (^j Vj +e{b-c) \yr]jf + e^{-c)b (V^rjj : 
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+ ^dt (1 + er;) Vf+e{d-a + der,) {VVj : VVj) + e^-a)d (V^V,- : 

< eCCabcdUj {U, ([/, + HUs + U!) + cj (t)(1 + C/«) {H^ + HU^ + U^)), (2.40) 

where {cj (t))^ is a sequence with {t))j G £’'(Z), having norm 1 and Cabcd is defined in p.39l) . 

Actually for the sake of simplicity, from now on we will not carry on the distinction between constants 
that depend on the a, b, c, d parameters and the ones depending on the dimension n and regularity 
index s. 

2.3 Another useful estimation 

At this point, working with the non-cavitation hypothesis: 

1 -I- erjo [x)>a> 0, (2-41) 

using estimate (|2.16l) and a bootstrap argument would be sufficient in order to obtain long time 
existence result similar to the one obtained in |12| (with some restriction on the value of e depending 
upon the initial data and a). However, proceeding in a slightly different manner we can avoid the use 
of (I2.41|) (although, all physically relevant data will verify it as 1 -|- er] represents the total height of 
the water over the flat bottom). We also stress out that the estimates established in this section are 
only available for the parameters a, b, c, d verifying m with the exception of the two cases (HH). 
Let us investigate the following quantity: 

hlLooC/^) =/l+/2, 

where 

r 2 h=eM\^^dtUf, 

\ 2h=sUfdt\\vh^. 

Owing to (12.111) we see that: 

(Ti +T2 + T3) + J divHj -t ce^ J ijV divVjVrjj^ 

< Ce^Us {U,Us + c, (t) + HUs + U^)) + Cs^ \\v\\l^ l|divkS||^. 

+ Ce^\\v\\l^ l|Vdivy,||^.||V7?,||^. 

< Ce^UjUs {UjUs + Cj (t) {H^ + HUs + Uf)) + Ce^U^U^ + Ce^U^U^ 

< CeUj {U, U^s + Cj (t) Us {H^ + HUs + U^)) ■ (2.42) 


Remark 2.1. Let us notice that the term c£^ f r/V div VjVrfj raises some important issues. In order 
to successfully estimate it (and thus in order to have the validity of (12.421) h we need the restriction 
dm) on the parameters a, b, c, d. The idea is that when c ^ 0, we must /lav^j 

sgn (b) + sgn (d) — sgn (a) > 2. 

In view of the fact that b + d > 0, it transpires that we must exclude the cases: 

a = (i = 0, b > 0 , c<0 and 

®One of the unknown functions y, V must have at least Sj^^-regularity level while the other one needs Sj^^-regularity 
level. 
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a = b = 0, d > 0, c < 0. 


Also, it is worth announcing that establishing (|2.42|) isn’t the only place where the restriction on the 
parameters is needed. As it will be soon revealed, in order to obtain local existence of solutions we will 
again have to bound —ce^f rjV divVjVrjj and the above considerations will have to apply. 

In order to handle I 2 we use the fact that the function t ^ 11^ (011^°° i® locally Lipschitz we get 
that a.e. dt ||?7 (i)llicxj exists and besides, a.e. in time we have 


\dt\\vmL^\ = 


lim 

s—yt 


bWllioo - his) 


t — S 


< lim 

S—¥t 


r]{t)-r]{s) 


t — s 


< \\dtri{t)\\^,^ < 119(7? (t)II a 


Thus, owing to (12.301) . (12.311) we get that: 


h < \eU]dt < CeU] {Us + HUs + U^) . 

Combining (12.421) and (j2.43[) we get that 

dt (I hILoo C/j) < CeUj {Uj {Us + HUs + U^) + (t) Us {H^ + HUs + U^)) . 

Finally let us observe that by adding (I2.16|) with (12.441) we obtain that: 

dtNf < CeUj {Uj {Us + HUs + U^) + c, (t) (1 + Us) {H^ + HUs + U^)) 


where Nj is the quantity defined in (I1.18|) . 


(2.43) 


(2.44) 


(2.45) 


2.4 The case 6 > 0, c < 0 and a = d = 0 


As pointed out in Section [221 , we are not able to establish (12.161) for the case (j2.15l) . However, if we 
proceed slightly different we can repair this inconvenience. Let us give some details of this aspect. As 
we have seen, the problem comes when estimating [ 7 ?, Aj] V div V (see (j2.231) 1. In order to bypass this 
problem, let us rewrite equation (12.61) in the following manner 


where 


-^dt (^J +eb\/^r]j : - ce J VdivF,V?7j 

—ce^ J Aj ( 77 V div V) Vrjj = T 5 

Ts = ce^ j {V^PjW : Vr]j) + ce^ J {V^rjjV : Vr]j) 

+ce^ ( / div V7?V7?j + / VVjVpVrjj) + ce^ / Aj {VWVr])Vrij — ce^ / RjfS/rjj 


(2.46) 


and 


Rjs = [W, Aj] + [y, Aj] V'^rj + [Vt?, A^] div V + [V 77 , A^-] VF. 
We add (I2.46|) with (12.31) and (12.121) in order to obtain 

^dt (^J Tyf + e{b - c) |V7?jf + e^b (-c) : V^tJj + (1 + et?) 
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I / A, (,V div V) V„ = £ / VTOr. + T. + T, + T.. 


Let us write that: 

„^2 / „v7 A . 


Thus, we get: 


'jVj = -ce^ 

J WrjArjjVj — cs^ J v^Vj div V,- 


= —ce^ 

J \7r]Ar]jVj — ce^ J Arjj [ry, Aj] div V — ce^ J 

ArjjAj (rydivT^) 

= —C6^ 

J \7r]Ar]jVj — cs^ J Arjj [ry, A^] div V + ce^ J 

\7r]jAj (Vrydivl/) 


VrjjAj (ryV div V) 


= —ce^ 

J VrjArjjVj — ce^ J Arjj [ry, A^] div V + cs^ J 

VrjjVrj div Vj 


Vrjj [Aj, Vry] div V + ce^ J VryAj (ryV div V) 


[ rjVArjjVj 

— ce^ / VryAj (ryV div V) 



= — ce^ J S/rjArjjVj — ce^ J Arjj [r], Aj] div V + ce^ J Vry div Vj 
+ ce^ J \7r]j [Aj, V 77 ] div V 

= —ce^ J WrjArjjVj — ce^ J Arjj [ry, A^-] div V — ce^ J V^rjjVrjVj 
— ce^ j S/^riVrjjVj + J yr]j[Aj,Vr]] div]/ 

c£^ (l|V?7||ioo \\Ar]j\\^2 \\Vj\\^2 + WArjjW^^Cj (t) ||V?7||j 


^ -OC I II VZ/lIrno Il^//Tllr9 II^Tllr^ T ||^//•^ll^9t.^ || V//|| || V \\ 

+ w^vh- I|v\-|L. \m\^, + \\v^vh^ w^vjWl^ mh. + 

+ l|Vry,|l^.c,(t)||V^ry||^._, II^IIb|,J 

< -cCs {Ups + Cj (t) Up) . (2.47) 

A similar reasoning as in (12.2311 together with estimates (12.221) . (12.291) . (|2.38|) and (12.471) gives us : 

\dt (^j ryf + e (6 - c) iVry^l^ + e^{-c)b (V^ry^ : V^ry^) + (1 + e?y) vf^ 

< CeUj {Uj {Us + HUs + Ul) + + Us) {H^ + UsH + U^ . (2.48) 

3 Existence and uniqueness 

We begin by defining what kind of solutions we are looking for: 

(I — ebA) dtT] + div V + as div At^ + sWVr] + e div (ryF) = 0, 

(/ - sdA) dtV + V?y + ceVAry + e^V \W\^ + sVWV + sWW + e^V |y|^ = 0, 

?y|t=o = Vo ) f^t=o = Vo¬ 


id 
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Definition 3.1. Let us consider a positive time T > 0 and ( 770 , Vq) G L^x and W G {H^ (R"))". 

A pair irj,V) G C([0,T],L^ x is called a solution to (11.101) on [0,T] if for any ((/>, 7 /’) G 

([0, T], 5 X 5”) and for all t G [0, T], the following identities hold true: 


I I I I 

J ivAI - £bA)dt(l))L2 + j (F, (/ + aeA)V0)i2 +ey {r],V {Wcf))^^2 + eJ {r]V,W(l>)^2 


= ivit), (I - sbA) (j){t ))^2 - (%, {I - sbA) (l)( 0 ))^2 


and 


J {V, (I - edA) dt-ip) 1^2 + j + ceA)div'ip)^2+eJ div!/;\ 

0 0 0 ' ' ^ 

e J (y, V*ll"-!/')i2 y (y,div(V'(8) W^))i2 +ey/^^,div?/;\ 

0 0 0 ^ 

= {Vit), (/ - shA) m)L^ - {Vo, {I - ebA) m)L^ ■ 


Let us state now the following local existence and uniqueness theorem which serves as an interme¬ 
diary result for Theorem [TJ 


Theorem 2. Let 0 , 6 , c,d be chosen as in m excluding the two cases dni), b + d > 0 , r G [0, 00 ] 
and s G R such that: 


s >—-|-1 or s =—-|-1 and r = 1 . 


(3.1) 


Furthermore, let us consider si, S 2 and S 3 defined by relation (11.121) and W G .Then, for all 

ido^Vo) G X with curlVo = 0, there exists a positive T > 0 and an unique solution 

(77, V)gC ([0, T], X if rex: or 

(77,F) G L- ([0,T],i?2% X n flC ([0,T],i?2V X (^2:^0;'')”) ^fr = 00, 

/3>0 


of equation (11.101) . Moreover, if we denote by T(? 7 o,Vo) the maximal time of existence then, if 
T {do, Vo) < 00 , we have that: 


lim 

Us (t) = 00 if r < 00 , 

(3.2) 

t^T{rio,Vo) 

lim sup 

Us {t) = 00 if r = 00 . 

(3.3) 


t->-T{rio,Vo) 


Of course, the local existence result of the above theorem is not optimal. For some particular 
choices of the a,b,c,d parameters we can solve the Cauchy problem for initial data in larger spaces, 
see for instance 0 , m, 0 . However, the lower bound on the time of existence is at most of order 
Initially, this is also the case of the solution constructed in Theorem [5] As mentioned 
above this time scale is not satisfactory from a practical point of view. However, as a by-product of 
the explosion criteria (|3.2I) - (I3.3I) of the solution of (11.101) and some refined energy estimate, we can 
improve the lower bound of the T ( 770 , Vq) thus establishing O (i)-long time existence. 

Having established all the estimates that we need, let us proceed by proving Theorem [21 
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Proof. We will use the so called Friedrichs method. For all m S N, let us consider Em the low frequency 
cut-off operator defined by: 

Em/ = (xB(O.m)/) ■ 

We define the space 

L‘L = [f ■■ Supp/ C B (0, m)| 

which, endowed with the H-H^ 2 -norm is a Banach space. Let us observe that due to Bernstein’s lemma, 
all Sobolev norms are equivalent on L^. For all m G N, we consider the following differential equation 
on 

( dt'n = Fm {v, V ), 

I dV = Gm{v,V), (3.4) 

[_ V\t=0 ~ Em?70) L|t=0 ~ EmVo, 

where (Fm, Gm) ■ ^ (^m)" ^ (Fm)" are defined by: 


Fm (? 7 , V) = —Em (^{I — ebA) ^ [(/ -I- aeA) div V + eWVr] + e div 
Gm iv. 1^) = -Em ((/ - edA)-^ [(/ + ceA) Vry + |v |F + W|"]) . 


(3.5) 

(3.6) 


It transpires that due to the equivalence of the Sobolev norm, (Fm,Gm) is continuous and locally 
Lipschitz on F^ x (F^) • Thus, the classical Picard theorem ensures that there exists a nonnegative 
time Tm > 0 and a unique solution (ry™,!/™) : ([0,rm],Fm x (Fm)")- Moreover, if Tm is hnite, 

then: 

^hm 11(77™, F-) 11^2 =00. 

Another important aspect is that because of the property = Em we get that the estimate obtained 
in (I2.11|) still holds true for {rjm, Vm), namely 

\9t (/ {vTT +e{h-c) |V77™|" + e\-c)h : V^ry™)^ (3.7) 

+ ^dt (^J |F™|^ +e{d-a) (VF™ : VF™) + e^{-a)d (VV™ : W™)^ 

+e J 77 ™ 77 ™ div F™ - J div F^Vry™ = Ti Fz -f T 3 


with Ti,T 2 ,T 3 dehned as in relations (12.4L (12.711 . (j2.10|l but with (ry™,F™) instead of (ry, F). Consid¬ 
ering C™ and up the quantities dehned in (iLTHll and (ITTTII with (ry™,F™) instead of (ry,F). Also, 
for the Tfs we dispose of the estimates (12.221) . (12.241) and (12.271) with F™ and F™ instead of Uj and 
Us and thus we gather that: 

Ti -H F 2 + T 3 < eCF™ (f™F™ + cj (t) -H HUP + {Upf'^ ) . 

Also, let us notice tha® 

-e J rj'^rip div F™ < e^G{upf Up. 

^^Observe that due to <|1.12|l at least one of rj, V has regularity level and thus, eventually by an integration by 

parts we can obtain the announced estimate. 
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Next, observe thalF^: 

J div VpVrif <eiC {UJ^f U^. 

We thus gather that: 

\dt (c/fUr + cj it) + iJ[/™ + 

and by a Gronwall-type argument we obtain that for all t G [0, T^]: 

t 

(t) < (0) +£^cj (upur + Cj (r) (h^ + HU^ + dr, 

0 


thus multiplying with 2-^® and performing an t'’’(Z)-summation, owing to Minkowsky’s theorem, we get 
that: 

t 

UT {t) < UT (0) + Ce^t + e^Cmax (1, H) J (u^ + dr. 

0 

We denote by H* = max (1, H). Thus, Gronwall’s lemma along with the explosion criterion for ODE’s 
gives us that: 

ln(l + ^j7^) ^ 

— TT-tr - ^ 

e 2 CH* 

where is the maximal time of existence for (1311). Also, due to the nature of the cut-off operator 
Em, it transpires that: 

ur (0) < Us (0), 


and consequently that: 


In (^1 -I- 

e^CH* 


< T* 


(3.8) 


In particular, for each time T > 0 such that the above inequality is strictly satisfied we obtain due to 
Gronwall’s lemma that for all t G [0, T] 


UT it) < 


j)-e^CH*T _ ^ 


(3.9) 


and thus, the solution {r}m,Vm) is uniformly bounded on [0,T]. Next, from relations (13.41) and (|3.5I) 
we get that idtv"^)mejq is uniformly bounded on [0,T] in Indeed, the first term is 


(/ - beA)-^ (J -k aeA) div G 


witfP^ 

S 4 = s -I- 2sgn(b) + sgn (d) + sgn{a) — 1 > s — 1. 

^^Again, because of 111. 1211 if c 7 ^ 0 then at least one of rj, V has regularity level while the other one has regularity 

level see Remark dril l. Thus, eventually by an integration by parts, we get the announced result. 

^^remember that b d > 0 
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In view of (|3I1), B 2 is an algebra and thus the last two terms of (I3.5|l are also at least in 
Thus, in view of the uniform estimates (13.91) . we get that uniformly bounded on [0,T] in 

Next, considering for all p G N a smooth function (j)p such that: 

f Supp{(l3p) C B {0,p+l) 

\ (t)p = l on B {0,p) 

it follows, in view of Proposition 16.51 that for each p G N, the sequence uniformly 

equicontiniuous on [0,T] and that for all t G [0,T], the set {<ppr]'^ {t) : m G N} is relatively compact 
in B 2 ~^■ Thus, the Ascoli-Arzela Theorem combined with Cantor’s diagonal process provides us a 
subsequence of F%nd a tempered distribution 77 G C ([0, T], S') such that for all (p G V (K”): 

inC([0,T],i?2V)- 

Moreover, owing to Proposition 16.31 and (13.91) we get that rj G (^ 2 ^) using interpolation we 
get that 

^ 077 in C ([0, T ], 

for all 7 > 0. Of course, by the same argument one can get V G ((^ 2 ^r)”) such that for all 
0G(T>(M”))”: 

inC([ 0 ,T],(B 2 V)”) 

for all 7 > 0. We claim that the properties enlisted above permit us to pass to the limit when 
777 — >■ 00 in the equation verified rf'' and P™. By the Patou property of Besov spaces we get that 
( 77 , P) G L“(i? 2 V ^ i^Tr) )■ It remains to verify that ( 77 , P) has the announced regularity. Suppose 
that r < 00. From 77 ’s equation we see that dtp G and thus, 77 G C ([0,T],which also 

implies that Sjp G C ([0, T], for fol j ^ The conclusion follows as the sequence of 
functions tends uniformly to 77 . Indeed, for all i,j G Z, we have 

Af(77 - Sjp) = 0 if ^ < j - 2 , 

therefore we have 


and as a consequence of 


which implies that p G C ([0, T], ^ similar argument shows that V G C ([0,T], When 

r = c», we know that for every positive /3, ^ continuously embedded in 00 ^rid repeating 

the above argument permits us to conclude that ( 77 , P) has the desired regularity. This completes the 

denoted the sake of simplicity. 
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proof of existence. Uniqueness, is a consequence of the following stability estimate: let us consider two 
solutions of (I1.10|) . (ry^,U^), and observe that the difference 

{Sv,6V) = {v^-7j\V^-V^) 

satisfies the following system: 

( (/ — ebA) dtbrj + div 5V + ae div A8V + eWV5r] + e div + e div {rf5V) = 0 

(/ - edA) dtSV + V6r] + ceVASri + eVW{SV) + eV{5V)W + £V{5V)V^ + eVV^SV = 0 
[ = 0 , 6V\t=o = 0 . 

(3.10) 

We consider U} = Ug and Ug = Ug see (11.171) . For the sake of simplicity we will 

prove stability estimates in the classical Sobolev space X = x with: 

f ri = sgn (&) - sgn (c), 

\ r 2 = sgn (d) - sgn (a). 

We endow X with the norm: 

Ih: ^ 11 ^ = \\v\\l 2 +£{b-c) WVrjWl^ - eHc || 

+ ll^lli^ +e{d-a) llVUIli, - e^da || 

Observe that due to the fact that si, S 2 are chosen so as to satisfy H1.12|) with s chosen as in (13.IL we 
have that x is continuously embedded in x (B^^) ■ Multiplying the first equation in 

(13.101) with g, the second equation with V and adding up the results we get that: 

{\\Sg\\l 2 +£{b-c) ||V( 5 ? 7||^2 + \\ 8 V \\\2 +£{d-a) ||V( 5 U||^ 2 ) + j ^dW 6 V 6 g 

+£c j AXSgSV < C^e {H + + U^) \\Sg, 5U||^ . 

Now, let us multiply the first equation in (13.101) with ceAg, the second equation with aeAV and adding 
up the results we obtain: 

]^dt (-c£ \\V 5 g\\l 2 - e^bc- ea || V< 5 U| 1^2 - e^da || 

+ eay X6gA6V + ec y div6VA6g 
<CV^{H + U^ + U^)\\6g,SVf^. 

Thus, adding up the above relations gives us: 

dt \\ 6 v, <CV^{H + Ul + U^g) \\Sg, <5U||^ . 

Hence, Gronwall’s lemma ensures the desired result. 

Proving the blow-up criteria is classic. Let us suppose that T (?7o, Vq) < c» and that 

limsup Ug (t) < 00 . 
t^T(-qo,Vo) 

Then, Ug (t) remains bounded on [0, T ( 770 , Vq)) say: 

Ug (t) < M, 
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for all t G [0,T {r]o,Vo)). We see that for any to G [0, T ( 770 , Vb)) we can construct, using the same 
method as before a solution to ( 11 . 101 ) with initial data (77 (to) ,V (to)) on a time interval that according 
to p. 8 l) satisfies the following lower bound: 

j^new _ ^ + ^ ln(l+^) 

° “ e^CH* ~ e^CH* 

Of course, choosing tg close enough to T ( 770 , Vg) we can obtain T”®*" > T ( 770 , Vg) such that gluing 
together the new solution with the ^)|[o to] in view of the uniqueness we get a contradiction on 
the maximality of T ( 770 , Vg). This concludes the proof of the announced result. □ 


4 The proof of Theorem [T] 


We are now in the position of establishing the announced long time existence result. Actually, we 
prove long time existence and uniqueness in the framework of the more general Besov space. More 
precisely, Theorem [T] is just a particular case of the following: 


Theorem 3. Let a,b,c,d as in OH) excluding the two cases (di), b + d > 0 . Let us take r G [1, c»] 
and s such that 


s> — + 1 ors= — + 1 and r = 1, 


with 77, > 1. Let us also consider si, S2 and S3 defined by (11.121) and W G . Then, we can 

establish long time existence and uniqueness of solutions (see DeRnition \1.2i) for the equation (jl.lOl) 
in i? 2 V ^ (-® 2 ^r) • Moreover, if we denote by T {rio,Vo), the maximal time of existence then there 
exists some T G [0, T (rjo, Vg)) which is bounded from below by an O (-)-order quantity and a function 
G : R —> R such that for all t G [0, T] we have: 


Us{v,v)<G{Us{m.Vo)) 
where Us{rj,V) is defined in relation (11.171) . 

Let us consider the unique maximal solution ( 77 , V) of (|1.10l) that we constructed in the proof of 
Theorem [5] Of course, the estimate (12.451) holds true for this solution thus we have that: 

dtNf < CeUj {U, {Us + HUs + U^) + cfit){l + Us) {H^ + UsH + U^)). (4.1) 

We recall that Nj is the following quantity: 


W = y (1 + e ll^ll +e{b-c){l + e II 77 II^„„) f 

+ j e^(-c)&(l+e||7?|licx,) (v277j 

+ J (l + e77 + e||77||^oo)F/+e(d-a + (ie77 + fie||77||^„„)(VV,- : VV,-) 

+ J e^{-a)d{l+eM\^^){y%:y%). 

and that we have: 

U, (t) < N, (t) <{l + 2e II 77 U, (t) (4.2) 
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thus we immediately get that 

dtN] < CeN^ {U, {U, + HU, + U^) + + U,) {H^ + U^H + U^)) 

and by time integration we get that: 

Uj (t) < Nj (t) < Nj (0) + Ce [ {Uj {Us + HU, + U"^,) + Cj(^)(l + U,) {H^ + UsH + U^)) dr. 

JO 

Multiplying the last inequality with 2^^ and performing a £’’(Z) summation yields: 

Us (t) <NQ + Ce [ {H^ + Us{H + H^) + U^ {1 + H) + u!) dr, 

Jo 

<No + CetH^ +eC{l + H + H‘^) [ {U, + U^) dr, 

Jo 

<No + CetH^ +eC{l + H^) [ {U, + U^) dr, 

Jo 


(4.3) 


for all t £ [0, T (r?o, Vq)) where 


No = 


(»)),« 


< (1 + 2e ||??o|Il<x.)^ ||(??o, %)\\b: 


Having established (14.31) we are in the position of finding a lower bound of order O (i) for the time of 
existence. Indeed by Gronwall’s lemma, and taking in account the explosion criterion of Theorem [5] 
we get that, supposing T ( 770 , Vb) is finite: 


limsup Us (t) = + 00 . (4.4) 

t^T{r]o,Va) 


Hence, via Gronwall’s lemma, we have: 




(iVo + £OT2r ( 770 , 1 ^ 0 )) 


<eC(l + i/2)T(77o,Ho). 


(4.5) 


If 

eCH^T{r,o,Vo) > No i.e. T ( 770 , Ho) > 

then, we have nothing else to prove. If this is not the case, the LHS member of (14.5|) is larger then 
i In and we get that: 


T{vo,Vo)> 



eC{l + H^) 


Thus, T (770, Vb) is always bounded from bellow by a quantity of order O (i). 

We now prove that on a O (i)-order time interval we dispose of uniform bounds for the solution 
of (I1.I0|) . Let us consider 

T* = sup {T G [0, T (770, Vb)) : Vt G [0, T ], U, (t) <2 {No + sCtH^)} . 
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Then, from (j4.3L we deduce that for all t <T* we have 

Us (t) < No + eCT^H'^ + eC {l + H^) (l + 4 (TVq + Us (r) dr 

and according to Gronwall’s lemma we get that: 

Us {t) < {No + eCT*H^) exp (^eT*C (l + H^) (l + 4 {Nq + eCT^H^f^ ) . 

Now, if there exists a /3 G (0, 2) such that 

exp (er*C (1 + H^) (l + 4 {No + eCT*H^y^^ < /3, 

then a continuity argument will lead us to the conclusion that T* = T {rjo, Vq) which will imply that 
T {vo, Vb) = 00 . Thus, for all t G [O, we get that 


Assume that 
Then, if 


Us {t) < ANo < (1 + 2eUs ( 770 , Vo))^ U, {r,o, Fq) • 
eT*C{l + H'^) (l+4{NQ+sCT*H^y'^ > ln2. 


T* > 


No 


eCH'^' 

we get that (14.6p holds true for all t G [O, ■ If (14.8p doesn’t hold i.e. 

No 


T* < 


then, combining (14.711 with (14.91) . we get that: 


T* > 


ecm 


In 2 


eC(l + i72)(i + i6Ar2) 

and thus, using the definition of T*, we get that for all t G 

/ 7/2 ln2 

Us (t) < 2 A^o + 


0 2 

’ eC{l+H^){l + 16Ng) 


(1 + i72) (1 + 16A^2) 


ln2 \ 

- r°+(l + 16iV2)) 

< 2 ("(1 + 2eUs (770, Vo)y Us (770, + 


In 2 


(I + I6C/2 (770,^0)) 


Thus, one can choose 


(4.6) 

(4.7) 

(4.8) 

(4.9) 


(4.10) 

(4.11) 


— mir, / (l+2|a:|) 2 X _ln2_ 

f (X) — mm < , c{i+h‘-‘)(i+iqx^) 

G{x) =max|2(^(l + 2|a;|)=a;+, 4 (1 + 2 |a:|) ^ a;| 

in order to conclude that T* > F (C/q) and for all t G [0, T*] we have: 

Us{t)<G{Uo). 

Remark 4.1. Theorem\^is just the restatement of the above result when r = 2. 

Remark 4.2. Observe that the above arguments allow us to derive similar uniform bounds for the 
quantities {2^^Nj^ 
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5 Final remarks 

The results of this paper generalize a part of the long time existence results for the abed Boussinesq 
systems that can be found in m- Theorem [3] has also a homogeneous counterpart i.e. one can replace 

the nonhomogeneous Besov space x with (^B 2 i ^ ^ (^ 2^1 ^ ^ definition 

and some basic properties about homogeneous Besov spaces see |31, Chapter 2. To our knowledge, 
B 21 n - 62 V) ^ (-^ 2^1 C i® the largest space in which one can prove long time existence for the 

abed systems. Let us also point out that as opposed to classical Sobolev spaces, working with Besov 
spaces enables us to attain the critical regularity index s = § + 1 - 

5.1 The remaining cases 

As we have seen earlier our method requires the restrictions (11.81) imposed on the parameters a, b, c, d 
in order to obtain the local existence theory as well as estimate (12.421) . see Remark l2.ll However, the 
general estimations (j2.16p (see Section ITM are valid for all the parameters (jl.4l) with b + d> 0. Using 
(j2.16l) and supposing that 

1 + e? 7 o > a > 0 , (5.1) 

we can obtain long time existence for any solutions of ( 11 . 101 ) in the cases: 

a = d = 0 , & > 0 , c <0 and 
a = b = 0, d > 0, c < 0 . 

However the Friedrichs method is not well-suited in this case because when establishing (I2.16p we 
multiplied the frequency localized equation of Vj with rjVj. Thus we would run into trouble when 
establishing the (j2.16l> -tvDe estimate for the approximations ( 77 ™, U™) introduced in Theorem [2l (13.41) . 
Nevertheless, one can imagine another strategy in order to bypass this inconvenience namely establish¬ 
ing a convergence scheme for (11.101) . In order to keep a certain homogeneity in this paper we will just 
give some brief details of how this strategy would work. Let us take IT = 0 to fix the ideas. Consider 
the linear systems: 

(/ - ebA) + div -h ae div -F e div = 0, 

(/ - edA) dtV'^+^ -)- Vry^+i -h ceVAjy^+i -h = 0, {C^+i) 

^m+1 o ^ T/-m+l o T/ 

V\t—Q — '^m+ 1^0 ? 

with Sm = Aj, and ( 7 y°,U°) = (0,0). We claim that it is possible to establish the following 

_ 7 <m —1 

(12.161) -type estimate: 

^dt (^J + e{b-e)\V7j^+^\^ +e\-c)bV^r]"^+^ + 

]^dt (1 -h £7?™) +e{d-a + der,”^) \VV"^+^\^ + s^{-a)d\7^V^+^ : 

<sPi[\\{rj^+\V^+^)\\g^ ,yi||V(77"^+i,U™+i)||3, ,e||V 2 (77™+i,U™+i)||^, )x 

P 2 ^ 11^ ^ £ IIV" ^""^IIbi, J ’ 

where Pi and P 2 are two polynomials of degree 2 with coefficients depending only on the a, b, c, d 
parameters, on the dimension n and on the regularity index s but not on e. Afterwards, Gronwall’s 
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lemma combined with some stability estimates and a continuity argument would allow us to obtain 
long time existence for (|1.10l) in the remaining two cases. Of course, we would use in a decisive manner 
the non-cavitation condition (15.11) . 

5.2 The a6c(i-system with smooth general topography 

As we have mentioned in the introduction abed-type models have been established in for a bottom 
given by the surface: 


{(x, y,z) : z = -h + eS{x, y)}, 

where S is smooth enough. In this case the abed-system, reads 

f {I — sbA) dtt] -j- div V ae div AV e div {{rj — S)V) = 0, 

I iI-edA)dtV-^yr]-\-eeyAr] + e^V\Vf =0, (5.2) 

[ ??|t=o = V\t=o = Vq. 

We claim that our method applies to this system with some minor modifications. Let us give a few 
details. For the sake of simplicity we will suppose that curl Vq = 0. By localizing the above equation in 
Fourier space and proceeding in the same manner as in Section [2] we can see that the terms preventing 
us from establishing long time existence are 

e J {y — S) r]j div Vj and — ee^ J {rj — S)\7 div Vj\7r]j. 

In order to repair this inconvenience we proceed as in Section [m the only difference being that we 
must multiply the localized equation of V with (rj — S) Vj rather than simply yVj. Thus, we can obtain 
an estimate similar to 112.16|) for the following quantity 

\dt (yj Vj + £ (6 - c) \yyjf -t- £^{-e)b : V^yj)^ 

+ ^dt (^J (1 + e{y - Sj) y/ + £ (d - a + de{y - S)) (VVj : VVj) + e^{-a)d {V^Vj : < eG{Us) 

where G should be some polynomial function with its coefficients not depending on e. Afterwards 
proceeding like in the rest of the paper, and taking S' G ^ with some s large enough, one can obtain 
long time existence results for (15.21) similar to those obtained in Theorem [TJ 
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6 Appendix: Littlewood-Paley theory 

We present here a few results of Fourier analysis used through the text. The full proofs along with 
other complementary results can be found in [3]. In the following if 17 C K" is a domain then V (17) 
will denote the set of smooth functions on 17 with compact support and S will denote the Schwartz 
class of functions defined on R". Also, we consider S' is the set of tempered distributions on R". 

Proposition 6.1. Let C be the annulus S R" : 3/4 < |^| < 8/3}. There exist two radial functions 


X S 27(73(0,4/3)) and tp G 27(C) valued in the interval [0,1] and such that: 

VCgK", x(?) + E‘^(2-^'^) = 1, (6.1) 

j>0 

VeGR”\{0}, = (6.2) 

2 < |j - j'\ ^ Supp((/5(2“'’.)) n Supp((/5(2“'^ •)) = 0 (6.3) 

/ > 1 ^ Supp(x) n Supp(v3(2“.’.)) = 0 (6.4) 

the set C = B(0, 2/3) + C is an annulus and we have 

|j-/|>5^2^Cn2^'C = 0. (6.5) 

Also the following inequalities hold true: 

V^GR", i<x2(e)+E‘^'(2”'0<l> (6.6) 

j>0 

V^gR-UO), i<^¥>^(2-^C)<l. (6.7) 

jGZ 


From now on we fix two functions y and ip satisfying the assertions of the above proposition. The 
following two lemmas represent some basic properties of the dyadic operators. 

Lemma 6.1. For any u G S' we have: 

u = Ajit in S' (R"). 

J6Z 

Lemma 6.2. For any u G S' and v G S' we have that: 

AjAiu = 0 if \j - 2| > 2. 


6.1 Besov spaces 

The following propositions are a list of important basic properties of Besov spaces that are used in the 
paper. 

Proposition 6.2. A tempered distribution u G S' belongs to i?|,, if and only if there exists a sequence 
(cj)- such that {2^‘^Cj)^ G F'(1) with norm 1 and an universal constant C > 0 such that for any j G Z 
we have 

IIAjmII^ 2 ^ C'cj. 

Proposition 6.3. Let s,s G R and r,r G [l,oo]. 
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• is a Banach space which is continuously embedded in S'. 

• The inclusion B 2 ^ ‘Z B^ f is continuous whenever s < s or s = s and f > r . 

• We have the following continuous inclusion Bf ^ C C rl^ fc L°°). 

• (Fatou property) If is a bounded sequence o/i3|^ which tends to u in S' then u G i?| ^ 

and 

< Iminf . 

• If r < 00 then 

lim IIM — <S'ju|| = 0. 

j—>-00 ^2,r 

Remark 6.1. Taking advantage of the Fourier-Plancherel theorem and using fF7\ one sees that the 
classical Sobolev space coincides with -B| 2 ' 

Proposition 6.4. Let us consider m G M. and a smooth function f : K." — >■ R such that for all 
multi-index a, there exists a constant Ca such that: 

veeM” \dViO\<Ca.ii + \^\r~'"'. 


Then the operator f (D) is continuous from i?| ^ to i?2 r™ ■ 

Proposition 6.5. Let 1 < r < 00 , s G R. and e > 0. For all 4> G S, the map u ^ (j)u is compact from 

Blr to B^ 27 - 

Proposition 6.6. Let s > 0 and 1 < r < 00 . Then _B| ^ n L°° is an algebra. Moreover, we have: 


|uu|| 




< 


INL~IHI« +M\b, IIHI, 


In particular, if s > ^ or s = ^ and r = 1, B 2 j. is an algebra. 

We end this section with the following result concerning a commutator-type estimate. For a more 
general form of this result and its proof see [3] page 116, Lemma 2.100. 

Proposition 6.7. Let us consider s G M, r G [1, 00 ] such that s > 1 + ^ or s = 1 + ^ and r = 1. Let 
{u,v) G X denote by 


Rj = [Aj,u]d°‘v = Aj {ud°‘v) — uAjd°‘v, 
where a is any multi-index with |q!| = 1. Then, the following estimate holds true: 


\Bi 


is the space of continuous bounded functions which decay at infinity. 
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